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Abstract 

In a central potential the usual resolution of the Schrodinger equation 
in spherical coordinates, whose solutions are of the form VE' = R(r) Y^ (9, tp) 
= [w(r)/r] Y^ (6, ip), consists in substituting in \l/ the solutions R(r) or 
u(r) of the radial equations considered as the radial parts of the Schrodinger 
equation. However, the solutions must be supplemented with the bound- 
ary condition u(0) = in order to rule out singular solutions. There is 
l/-\ still no consensus to justify this condition, with good reason. It is based 

f*"*; _ on a misunderstanding that comes from the fact that the radial equation 

in terms of R(r) is derived from the Schrodinger equation, and the radial 
equation in terms of u(r) from the former, by taking the Laplacians of *!/ 
and of R(r) in the sense of the functions. By taking these Laplacians in 
the sense of the distributions, as it is required, we show that the radial 
equations are derived from the Schrodinger equation when their solutions 
are regular, but not when they are singular, so that the equations need 
not be supplemented with any supplementary condition such as u(0) = 0. 



1 Introduction : The usual resolution of the SE 
in spherical coordinates. 

In a central potential V(r), to which we shall confine ourselves, and which will 
be denoted by V r , the operators H, L 2 and L z all commute with each other, so 
that the eigenfunctions of H, or the solutions of the Schrodinger equation (SE) , 
are of the separable form 

tt(r, 9,<p) = R (r) Yf (9, V ) = H^ r /x ^ ^ (1) 
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where the factor 1/r is introduced for convenience, and the spherical har- 
monics iy (0, ip) are eigenfunctions of the square of the angular momentum 



L 2 y/(0,p)=*(*+l)ft 2 y/ , (0, ¥ >) 

The usual resolution of the SE in spherical coordinates, 



(2) 



H^ 



h 2 1 d 2 



2m r dr 2 2mr 2 



V r 



R(r)Y?($,<p) = ER{r)Y?{$,<p) (3) 



is based on two assumptions. 

The first assumption is that in any case the radial and angular parts of the 
Laplace operator separately act, and hence that one may substitute (2) in the 
expression of the Laplacian of \f which then amounts to a radial part times 
Yf {6, ip), that is, 



A* = 



lcP_ L 2 
r dr 2 h 2 r 2 



R(r)Y?{e,<p) 



1 d 2 £(£ + 1) 
r dr 2 r 2 



R{r)Y?{6,y) 

(4) 



as also the Hamiltonian of "f, 

h 2 i d 2 



Ht> = 



2m r dr 2 



2mr 2 



h 2 1 d 2 £(£ + i)h 2 

2m r dr 2 2mr 2 



V r 



Vr 



R{r)Y?{0,<p) 



(5) 



The SE then reads 



ff* 



2m r dr 2 



£{£ + l)h 2 
2mr 2 



Vr 



R(r)Y?(6,<p)=ER(r)Y?{0,<p) 



which gives the radial equation 



h 2 1 d 2 £(£+l)h 2 

2m r dr 2 2mr 2 



Vr 



R(r) = ER(r) 



(6) 



(7) 



The second assumption is that by substituting u(r) and multiplying both 
sides of (7) by r, the latter amounts to the radial equation 



2m dr 2 



£(£+l)h 2 
2mr 2 



V r 



u(r) = E u (r) 



(8) 



The resolution of the SE (3) in the three variables r,0,ip, is thus replaced 
by the resolution of (7) or (8) , considered as the radial parts of the SE, or the 
radial SE, depending only on the variable r. The solutions, substituted in \f, 
are written in the form of series, R (r)= r s T,k>oOkr k , and u (r)=r s+1 T,k>oakr k , 



by considering that for potentials which are not more singular at the origin than 
1/r, the indicial equation gives the two roots s = I and s = - (£+1). For t > 
only the solutions given by the first root are normalizable. In the radial case [i 
= 0) Yq° —1/\^4tt, the equation (8), which amounts to 



2m dr 2 



u (r) = Eu (r) (9) 



is called the "reduced Schrodinger equation". Its normalizable solutions are given 
by both roots, but for s = -1, as u(0) = a ^ 0, R(r) and 5" behave at the 
origin like 1/r, the solutions are ruled out as inadmissible, so that the equations 
are supplemented with the boundary condition 

u(0) = (10) 

The question has always been to justify this supplementary condition. If 
these functions \& are considered indeed as solutions of the SE , the question is 
to justify that normalizable solutions of the SE are not admissible, and if these 
functions are not considered as solutions of the SE, the question is to justify that 
u[r) is nevertheless a solution of the reduced SE. This is why various arguments, 
often conflicting, have been proposed to justify this condition, but no one is 
satisfactory. It comes from the fact that the usual assumptions, and hence the 
equations which are based on them, hold in R 3 /{0}, but not necessarily in R 3 . 

2 The Hamiltonian in R 3 . 

Let us briefly recall that if $ € C°°(R 3 ) and s is a negative real, \&= r s Q, which 
is of C°°(R 3 /{0}), but not of C°°(R 3 ), defines in R 3 a distribution called "pseud- 
ofunction"[l,2], denoted by Pf.^f, and its Laplacian in R 3 is the Laplacian in the 
sense of the distributions. In the case $= 1 the theory of distributions shows 
that depending on whether or not s is an odd negative integer, its Laplacian 
has a different expressionjl], which can be written[3] 

(4r>+ 1W 3 / 2 
APf.r s = 8(8 + l)Pf-r s - 2 + x P C P APS C p = -- '' '' 



22p-i p i r( p + 3/2) 

fl if P^N s + 1 

Xp = \0 if P {N P= -— (11) 

where A p is the iterated Laplacian, and where the symbol Pf. is used be- 
cause depending on the value of s we have an ordinary function r s , or a pseud- 
ofunction Pf.r s . The significance of the term s(s+l)Pf.r s ~ 2 is the following. 
Just as r s , which is of C°°(R 3 /{0}), but not of C°°(R 3 ) if s < 0, defines in R 3 
the pseudofunction Pfr s , its Laplacian in R 3 /{0}, which is the Laplacian in 
the sense of the functions, Ar s = s(s + l)r s ~ 2 , defines in R 3 the pseudofunction 
Pf.Ar s = s(s+l)Pf.r s ~ 2 . Eq.(ll) therefore also reads 



APf.r 8 = Pf.Ar 3 + x P C P A p S 



(12) 



The symbol Pf. can be dropped in the case of pseudofunctions less singular 
at the origin than l/r 3 [l], but in the case of the operator Pf.A only if it is 
equivalent to A, e.g. when it acts on a function of C°°(R 3 ). Thus, Pf.(l/r) 
reads as 1/r, but P/.A(l/r) is different from A(l/r), for s = -1, (12) gives 



A(-) = Pf.A(-) + C S = -And 
r r 

The pseudofunction 

Pf.R (r) = Pf.r s Y^ a kr k 
where s is any real, is written 



Pf.A(-) = 
r 



(13) 



(14) 



fc>0 



Pf.R (r) = Pf.[r s S e (r)] + Pf\r s+1 S d (r)] 

where S e (r) = J2k>o a2kr 2k a,nd S d (r) = J2k>o a sk+i r 2k are of C°°(R 3 ), so 
that APf.r s S e (r) and APf.r s+1 Sd (r), which are the Laplacians of the product 
of a distribution by a function of C°°(R 3 ), are given by the usual expression of 
the Laplacian of the product of two functions. It enables to show that [3] 



APf.R (r) - Pf.AR (r) = Q s (5) 



(15) 



where 

Qs {6) = X P 



(C p + 2L p r—)S e (r) 



& p S + Xq 



(C q + 2L q r-)S d (r) 



A q S 



, C p -i C p a + 1 s + 2 

L v = — r-i- G_i =0 p = a = 

p 8p(2p+l) 4 1 y 2 y 2 

Consider now the pseudofunction 

P/.tf = P/.P (r) K/* (0, p)= P/.r s E fc >o «fcr fc 7/ (0, p) 
Substituting 

AF/.[r s ^ {0,<p) S e (r)] = AP/.{r s "'[/F/ (0 >¥ >)S e (r)]} 



(16) 



(17) 



AP/.[r s+1 r; (0, ^)S d (r)] = APf.{r s+1 ~Y Yf (9, <p)S d (r)]} 

AP/.\& is given by the sum of the two right-hand sides, where r £ T/ 1 (0, <p)S e (r) 
and r F/ 1 (0, <p)S d (r) are of C°°(R 3 ), so that these two Laplacians are the prod- 
ucts of a distribution by a function of C°°(R 3 ), given by the usual expression of 
the Laplacian of the product of two functions. It enables to show that [3] 



AP/.tf - P/.A* = Q 8i , (8) 



(18) 



where 



,/ (<*) = X P 



(C p + 2L p £ + 2L p r—)S e (r) 



d_ 
dr' 



r l Y£ (0, if) A p 5 



-Xq 



{C q +2L q £ + 2L q r—)S d (r) 



P 



s + 1- 



q = 



r £ Y^(e,if)A q S 

s + 2-£ 



(19) 



2 * 2 

These results hold in Cartesian and spherical coordinates. The Laplacian of 
Pf.R(r) therefore reads 



(l-f^r)Pf.R (r) = Pf.(L* r )R (r)) + Q 8 (6) 

r ar z r dr z 

or, in terms of Pf.u(r) = P/.r s+i £fc> afcr fe 

( 1 -^r)PfM=Pf.( 1 - d ^l ) + Qs (5) 
r ar z r r dr z 

where Q s (5) is given by (16). The Laplacian of Pf.il/ reads 



AP/.# = Pf. 



(i* _ffi + i) )J1(r) y r (,, rt 

r dr r 



Pf- 



\,d 2 £{£ + !). . ,.,„,. , 



+ Qs.i (5) 



and the Hamiltonian of $, 



(20) 



(21) 



(22) 



HPf.V = Pf. 



h 2 l d 2 £(£ + i)h 2 



2m r dr 2 



2mr 2 



+ V r )R(r)Yf(9,ip) 



-h^ 



Pf- 



l h 2 d 2 ] £{£ + \)h 2 



r 2m dr 2 2mr 2 

where Q s j (S) is given by (19). 



+ V r )u (r) Y? (0, p) 



2m 



Q,,e (5) (23) 



3 The usual assumptions in R 3 . 

In R 3 /{0} Q s (8) and Q s j (8) vanish (and there is no symbol Pf.), so that the 
equations (4) to (9) hold. Eqs.(20) to (23) show that a necessary condition 
for the usual assumptions to hold in R 3 is that Q s (5) = = Q s j (8) (it is 
not necessarily the case). In order to determine Q s (8), instead of (16), by 
substituting (12) and (14) we write APf.R (r)in the form 

APf.R(r) = J2*kAP.f.r k+s = ^> fc [p.f.Ar k+s + Xp{k) C p{k) A^8 

k + s + 1 



k>0 fc>0 



= Pf.AR(r) + J2^X P (k)C p{k) A^S 

fe>0 

which gives 



k + s + 1 



P 



Q s (8) = J2 a kX p(k)C p{k) A p{k) 8 

fc>0 

As p G TV, only if s is an integer < -1, and k < - s -1, the sum amounts to 
Q s (5) = J2 a kX P (k) C m A^S p=- *±1±1 (24) 

fe=0 

where all the terms of the same parity as s are zero. As C p ( k ) ^ =^ A P ^S 
(p(k) € N ), and as the A p ( fc ) (S) are linearly independent, Q s (8) ^ if there 
is one term of the sum at least which is non zero. Hence Q s (5) is non zero, 
and the second assumption on which the passage from the equations in terms 
of R(r) to the equations in terms of u(r) is based does not hold if there is at 
least one coefficient a k ^ for which k + s is an odd negative integer, that is, 




(,*)*() — =^/ __ 2p _ 1 P^ N ( 2 5) 

It is the case in particular if there are two non zero successive coefficients, 
a k =/= 7^ a/c+i, or/and if Pf.R(r) behaves at the origin like r s , where s is an 
odd negative integer (since a =£ 0). Thus, the passage from the equations in 
terms of R(r) to the equations in terms of w(r), in particular from (7) to (8), 
which is usually based on the relation 

,1 d 2 , . , ,ld 2 .u(r) ld 2 u(r) , s 

(~^r)R (r) = (-^r)— = --^5- (26) 

does not hold if the condition (25) is satisfied. It is the case for the solutions 
normalizable such that u(0) = a 7^ 0, and thus such that R(r) behaves at the 
origin like 1/r , since for s = -1, Q s (8) = a C 8, (21) gives 



Ad 2 , , . Ad 2 ,u(r) 1 d 2 u(r) , . „ , . 

{-—r)R (r) = (-— rii - — -^ - 4tt«(0 <J (27) 

r dr z r dr z r r dr z 

Note that this relation, very useful, can be obtained only with the help of 
the well-know relation A(l/r) = -AttS. If s € N indeed, except for s — 1, 
r s <G C°°(R 3 ), Ar s is the Laplacian in the sense of the functions, and as Ar = 
A [r 2 (1/r)] , the latter, which is the Laplacian of the product of a distribution 
by a function of C°°(R 3 ), is given by the usual expression of the Laplacian of 
the product of two functions, that is, Ar = 2/r, so that if s <G N, in any case 
Ar s = s(s + l)r s ~ 2 . Hence, if u(r) = Y> k > a k r k , 



7V ( 7* ) CI 

A— — = — + S fe >ia fe Ar fe " 1 = -Ai:a 5 + Z k>1 a k k(k - l)r fe " 3 
r r 

= -47ra (5 + S fe > afefc(fc - l)r fe - 3 = 11^1 _ 4^(0)5 
Similarly, the Laplacian of Pf.ty given by (17) reads 

AP/.tt = E fc >o «^P/.r fc + s y/ (0, p) 
By substituting 5 e (r) = 1 and Sd (r) = in (18) and (19), we have 

APf.r s Y? (9, ip) = Pf.Ar 8 Y$ (6, i p)+ Xp B^C p r e Yf (6, tp)&8 

s + l-l 

P= 2— 

where B Lp C p = C p + 2L p £ = [1- 21/ (ip +l)\C p . As I e N, B t , p = only if 
p i TV, so that XpB e , p C p ^ if p e A, and = if p <£ N. As r e F/ (0, p) is a 
homogeneous polynomial of degree £, r e Y\f (0, ip)A p (6) ^ only if 2p > I, and 
then only if s< -1, so that 

X P B Wp C p /Y^(e^)A^l^ if S ~lZ~T~] 8<-l,peN 

1 = if s — t^—2p—\ 

In particular for s = - (£+1), AP/.r-< l+1 Y>? (0, p) ^ 0. Hence 

AP/.* = ^ a fc [P.f.Ar fc +* y/ (0, <p) + x P (k)B^ p{k) C p{k) r e Y» (0, p)A^S] 
o 

= P/.A* + Q s ,i (5) P = 



k>0 

k+s+i-e 



where 



k+s+l-£ 



P 
with 



where s < — 1, p£ JV. As fc < - s - 1 (since 2p > £), the sum amounts to 

-s-l 

Q a>I (<J) = ]T a fc x pW ^( fc )C p(fc) ^y/'(0^)AP( fe )<5 

fe=0 

,..i±i+i^ (28) 

where the terms corresponding to even, resp. odd, values of k are zero , if 
s and i have the same parity, resp. different parities. In any case, Q S} i (5) ^ 
if one term of the sum at least is non zero. Hence Q s j (5) is non zero, and 
the first assumption on which the usual resolution of the SE is based, namely, 
the Laplace and Hamilton operators separate the variables so that AW and H& 
amount to a radial part times Yjf (6, ip) , does not hold if there is at least one 
coefficient a^ ^ for which k + s - I is an odd negative integer, that is, 



Q,AS)¥>0 <=► 3 C/° , , P^ N ( 29 ) 

I k + s — t = —2p — 1 

It is the case in particular if there are two non zero successive coefficients, 
a k 7^ 7^ a/s+i, or/and if Pf.^f behaves at the origin like r s Y£ (6,<p), where 
s - £ is an odd negative integer. Notably, if P/.W behaves at the origin like 
r -(i+ 1 ) y^ {9, ip), the first term is non zero, 

aoXiB^dr" Y£ (9, ip)A e 5 + 
and hence 



k—l 

-(i+i),i = J2 a ^P^ B ^(k)C P ( k )r e Y^(9,ip)A^S^0 p = ™^- (30) 

fc=0 



4 The radial equations and the radial Schrodinger 
equation. 

One usually considers that the solutions of (7) are given by the two roots of the 
indicial equation s = £ and s = - (£+ 1), and hence that they behave at the 
origin, either like r , or like P 'f.r~^ +1 > . As they are, either regular functions 
R{r), or singular functions, or pseudofunctions Pf.R(r), (7) should be written 



h 2 l d 2 



£(a + i)h 2 



V r 



2m r dr 2 2mr 2 

By substituting (20), the left-hand side reads 



Pf.R(r) = EPf.R{r) 



(31) 



Pf- 



2m r dr 2 
2m r dr 2 



£{£+l)h 2 
2mr 2 

£{£l + l)h 2 
2mr 2 



Vr 



Vr 



Pf.R{r) 



, s h 2 

y ' 2m 



.(*) 



(32) 



Eq.(31) can thus be satisfied only if the left-hand side involves no expression 
containing S, and then if Q s (S)= 0. Eq.(25) shows that if s = £, as k + £ ^ -2p 
-1 Vfc and Ml, Q e {S)= W. If s = - (£+ 1), Q-( t +i) ^ if there is at least a 
coefficient a^ ^ for which (£ — k)/2 G TV. It is the case if there are two non 
zero successive coefficients, ^ ^ O5/ a-k+i, or/and if £ is even . In this case 
(7) or (31) has no solutions which behave at the origin like Pf.r~( i+1 >). The 
solutions of (7) given by the root s = - (£+ 1) are solutions of (7) in R 3 {/0}, 
but not necessarily in R 3 . 



The radial equation which has solutions which behave at the origin like 



or Pf.\ 



-(M-l) 



, is the extension in R of (7), that is, the equation 



PI- 



2m r dr 2 



£(£ + l)h 2 
2mr 2 



V r 



R{r) = E Pf.R (r) 



(33) 



and similarly, the radial equation which has solutions which behave at the 
origin like r l+1 or Pf.r~ e , is the extension in R 3 of (8), that is, the equation 



PI- 



2m dr 2 



i)h 2 



2mr 2 



Vr 



u{r) = EPf.u(r) 



(34) 



(In other words, Pf.R(r) is a solution of (33), and Pf.u(r) a solution of (34), 
if in R 3 /{0} R(r) is a solution of (7), and u(r) a solution of (8), so that the 
solutions of (33) and (34) are given by the two roots s— £ and s = - (£+1)). By 
substituting (33) in (23) we obtain 



HPf.y = Pf. 



2m r dr 2 



£{£+l)h 2 
2mr 2 



Vr 



R{r)Y?{0,<p) 



EPf.[R(r)Yf(6,<p)}-—Q s ,e(6) 

2m 



—( 

2m 



(35) 



and by substituting (34) in (23), 



HPf.V = Pf. 



l ( _&_£_ 

r 2m dr 2 



e(£ + i)h 2 

2mr 2 



+ V r )u (r) F/ (6, <p) 



2m 



Q.,i (S) 



= EPf. 



u(r) 



r,"(o,v) 



2m 



:(S) 



(36) 



It shows that depending on whether or not Q s ^ (6) is zero, and hence de- 
pending on whether or not the condition (29) is satisfied, the radial equations 
are derived from two different equations. 

We can therefore conclude. When the solutions of the radial equations are 
given by the root s —£, they are regular functions, the symbol Pf. is useless, as 
k 7^ -2p -1 Vfc, Qt t t (5)= W, the radial equations are derived from the SE, the 
equations (4) to (9) hold in R 3 . 

When the solutions are given by the root s = - (^+1), they are singular func- 
tions, or pseudofunctions, and as Q-n+t) i{S) =£ 0, eq.(30), the radial equations 
are not derived from the SE, but from the equation 



HPf.ty = EPf.y ( 

2m 



-(e+i)j(o~) 



(37) 



In this case the radial equations are not the radial parts of the SE, this is why 
their solutions substituted in Pf.^f do not give solutions of the SE, but solutions 
of (37). For £ > these solutions are not normalizable, and then generally not 
taken into account, so that the fact that the radial equations (7) and (8), that 
is, (33) and (34), are not derived from the SE, and then do not hold, is of no 
consequence. 

The problem occurs in the radial case owing to the fact that the solutions 
are normalizable, and then must be taken into account. As the symbol Pf. can 
be dropped for \& and R(r) which behave at the origin like 1/r , but not for the 
operator Pf.A, (33) reads 



Pf- 



h 2 id 2 

2m r dr* 



R(r) = ER{r) 



(38) 



It is derived from the equation (35) for £= 0, that is, 



if* = 



1 



Pf- 



_]?_ld?_ 

2m r dr 2 



R(r) 



i(0)S 



ER{r) 



-u{0)8 



(39) 



As u (0) = a 7^ 0, u(r) is a regular function, Pf(d 2 u(r) /dr 2 )=(d 2 u(r) /dr 2 ) , 
(34) is written in the form (9). The latter is derived from (36) for £= 0, which 
reads, in agreement with (27), 



10 



i?* = 



1 1 


' h 2 d 2 

2m dr 2 


u(r 


y/4n r 




i E <r) n 2 ^ u 



-u(0)S 



4n 



u(0)S 



(40) 



Thus, usually (9) is considered as derived from the SE, and hence as the 
reduced SE, and the question is to justify that its solutions substituted in \P 
give solutions of the SE, and are admissible, when u(0) = 0, but not when u(0) 
7^ 0. Eq.(40) shows that in fact (9) is derived from the SE and is the reduced 
SE when u(0) = 0, but that when u(0) ^0 (9) is not the reduced SE, it is not 
derived indeed from the SE, but from the equation 



HV = EV 



-u(0)S 



(41) 



This is why the solutions of the radial equation (9) substituted in "J give 
solutions of the SE when u(0) = , but not when w(0) = a ^ 0. In this case 
they give solutions of (41). 

Eqs.(20) to (23) show that the usual misunderstandings occur when the 
Laplacian of W in (4) and (5), and hence in the SE (3), and the Laplacian of 
R(r) in (7), are not taken in the sense of the distributions, but of the functions. 
The radial equations (7) and (8) are then considered as derived from the SE 
regardless of their solutions. It can be summarized as follows : the solutions of 
the radial SE satisfy the condition u(0) = 0, they are normalizable, and they 
are solutions of the radial equation (9), but the normalizable solutions of (9) 
such that m(0) ^ are not solutions of the radial SE. This is why the boundary 
condition u(0) = is aimless. 



5 Concluding remarks 

The radial equations are usually written, either in terms of R{r), or in terms 
of u(r), the two equations being considered as equivalent. In fact it is not the 
case, and not only because the latter is derived from the former. Eq.(32) shows 
that when Q_^ +1 )(5) ^ 0, the solutions of (33) which behave at the origin like 
pj r -(t+i) are na (- so l u tions of (31), but of the equation 



fV_lf_, 

2m r dr 2 



£(e + i)h 2 

2mr 2 



Vr 



Pf.R{r) 



EPf.R{r) 



2m 



Q-(t+i)(6) 



(42) 



As the radial equations (33) and (34) are equivalent, to avoid any confusion 
between (7) or (31) and (33), it may be preferable to write the radial equations 
in terms of u (r) , all the more as in the case of the normalizable solutions these 
equations are written without the symbol Pf.. Besides, if u(r) is a function 
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singular at the origin, or a pseudofunction Pf.u(r), the term (d 2 /dr 2 )Pf.u(r) 
is not defined in the theory of distributions (and a fortiori of the functions), 
so that (9) holds in R 3 only for regular functions, which is just the case of 
normalizable solutions to which we confine ourselves. 

What we called the "usual resolution of the Schrodinger equation" is the 
process [4] in which the solutions of this equation are obtained by substituting 
in vf, either the solutions R(r) of (7), or the solutions u(r) of (8), where (7) 
has been derived from the SE (3), and (8) from (7), by taking the Laplacians 
of \I/ and of R(r) in the sense of the functions, while considering that there are 
singular solutions which must be ruled out. In other words, in no handbook 
of quantum mechanics, (nor article to my knowledge) the radial equation (7) 
is derived from the SE (3), and the radial equation (8) from (7), by taking the 
Laplacians of \P and of R(r) in the sense of the distributions, as it is required. 
Singular solutions are sometimes ruled out by considering that if R(r) is a 
singular solution of (7), "f cannot be a solution of the SE since A<3/ involves 
an expression containing 5. But if A\P involves an expression containing 6, it 
shows that (7) which does not involve itself any expression containing S, cannot 
be derived from the SE, in this case indeed R(r) is not a solution of (7), as it is 
assumed, but of (42). 

It is thus because the Laplacians are not taken in the sense of the distri- 
butions that the two usual assumptions on which this resolution is based, and 
then the equations that they imply, hold in R 3 /{0}, but not necessarily in R 3 . 
It is what constrains to impose the supplementary boundary condition u(0) = 
0. Justifying this condition has always been a requirement (See e.g. [5]), and 
remains a present-day problem (See e.g. [6]). The various and often conflicting 
arguments proposed to justify this condition are based on this misunderstand- 
ing, they never led to a consensus. Insofar as they are a part of the usual 
resolution of the SE, these arguments will be also revisited. 
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1 Erratum. In [3, p.17] we wrote , 2 L t l + C t = C £ _i/4(2£ +1) + C e /2, instead of C<>_i/4(2£ 

-l) + Ci{i-e/2). 
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